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Abstract: We apply the on-shell tree- level recursion relations of Britto, Cachazo, Feng 
and Witten to a variety of processes involving internal and external massive particles with 
spin. We show how to construct multi-vector boson currents where one or more off-shell 



X 



vector bosons couples to a quark pair and number of gluons. We give compact results 
for single vector boson currents with up to six partons and double vector boson currents 



with up to four partons for all helicity combinations. We also provide expressions for single 
vector boson currents with a quark pair and an arbitrary number of gluons for some specific 
helicity configurations. Finally, we show how to generalise the recursion relations to handle 
massive particles with spin on internal lines using gg — > ti as an example. 



1. Introduction 



Powerful new formalisms have been recently developed for the calculation of on-shell quan- 
tities such as scattering amplitudes in gauge theories. Most striking are the MHV rules of 
Ref. [1] and the Britto-Cachazo-Feng-Witten (BCFW) recursion relations of Refs. [2, 3]. 
Applications of these new formalisms together with the classic unitarity based approach of 
Refs. [4, 5] have led to a dramatic progress in calculations of amplitudes. 

The 'MHV rules', originally proposed in [1], led to simple and compact calculations of 
tree level amplitudes [6-12]. The method has also been generalised to include external Higgs 
bosons and massive vector bosons [13-15] as well for investigating collinear factorisation 
[16, 17]. 

The on-shell BCFW tree-level recursion relations were first formulated for massless 
particles [2, 3, 18-21]. The proof of the new recursion relations [3] relies on simple analytic 
properties of the amplitudes and have recently been generalised in three ways. First, the 
massless on-shell recursion relations in gauge theory were extended to include gravity [22, 
23]. Second, a new version of recursion relations was developed in Refs. [24, 25] to calculate 
all finite one-loop amplitudes in non-supersymmetric QCD. Ref. [26] further generalised 
this approach to compute the cut-nonconstructible parts of divergent one-loop gluonic 
amplitudes. In a third development, Ref. [27] generalised the BCFW recursion relations to 
include massive scalar particles at tree level. The goal of this paper is to further extend 
the recursion relations to include massive particles with spin such as vector bosons or 
massive quarks in a natural way. We will show that the recursion relations lead to new 
compact formulae for scattering amplitudes that are somewhat simpler than the previous 
calculations of Berends, Giele and Kuijf [28, 29]. 

Much progress has also been made at one loop level. The MHV rules have been suc- 
cessfully applied to compute supersymmetric amplitudes [30-36] as well as new calculations 
using improved unitarity methods [37-47]. The new formalism has been largely inspired 
by Witten's idea of a duality between supersymmetric Yang-Mills and a topological string 
theory on a twistor space [48]. For reviews of recent developments see [49, 50]. 

Our paper is organised as follows. In section 2 we review the recursive formulation of 
gauge theory amplitudes. Section 3 defines our model and reviews colour ordering of the 
amplitudes. We provide new compact results for amplitudes with up to six partons and 
a single vector boson in section 4 and with up to four partons and two vector bosons in 
section 5. New compact expressions for n-point NMHV currents with a single vector boson 
are also given in section 4. Finally we demonstrate how to use recursion relations including 
propagating massive particles with spin in section 6. Our findings are briefly summarized 
in section 7. An appendix stating our spinor conventions is also included. 
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2. On-Shell Recursion Relations 



Here we give a brief summary of the on-shell recursion relations [2] for tree-level amplitudes. 
These recursion relations follow from general quantum field theory arguments [3] and we will 
present them in the form [27] which naturally incorporates massless and massive particles. 
We consider the colour-ordered partial amplitudes A = A(pi, . . . ,p n ), in which the coloured 
particles come in a definite cyclic order 1,2, ... ,n. These amplitudes are obtained by 
stripping away the colour factors from the full amplitude, hence, they depend on the 
kinematic variables, momenta and helicities, pk and hk only. 

The formalism relies on choosing two particles in the amplitudes to be shifted by a 
complex vector to be specified below. We label these marked particles i and j. The on-shell 
recursion relation for a tree-level amplitude takes the form: 

Anipi, ■ ■ ■ ,Pn) = ^^(Pr, ...,Pi,...,p s , ~P h ) p2 ]_ m 2 



7 1 lltp 

partitions h r 

x A R (P~ h ,p a +i, . . . ,Pj, . . . ,p r -l) , (2.1) 

where P = p r + . . . + pi + . . . + p s and the "hatted" quantities are the shifted on-shell 
momenta. The summation in (2.1) is over all partitions of n external particles between the 
smaller amplitudes on the left, Al, and on the right, Ar, such that pi is on the left, and 
Pj is on the right, and also over the helicities, h, of the intermediate state. 

In this paper we will always choose both marked particles i and j to be massless. We 
shift two massless momenta 1 pi = and pj = of the marked particles by r\ = 
such that the shifted momenta are 

Pi=Pi + z\j)\i] , p1 = 0=p1, (2.2) 
ft = Pj - z\j)\i] , f J= Q=p), (2.3) 
P = P + z\j)\i], P 2 = m 2 P . (2.4) 

Here mp is the mass of the particle on the internal line. Equation (2.4) determines the 
variable z as 

P 2 — m 2 P 2 — m 2 
2P-r) (j\P\i] ■ { ■ ) 

For the particles i, j momentum shifts above are equivalent to shifting the spinors as follows: 

\i] = \i] , |?) = \i) + z\j) , (2.6) 
li) = |J> , \3] = \J] ~z\i]. (2.7) 



A direct proof of (2.1) was given in reference [3] based on the analyticity structure of 
the meromorphic tree- level amplitude A n (z) in the complex z-plane. For equation (2.1) to 

1 Hclicity spinors \i] are denned in the Appendix along with our conventions for forming the inner 
product. 
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hold it is essential that A n (z) should have no poles at infinity. In general it is known that 
whether or not these "boundary" contributions are present (i.e. whether or not (2.1) is 
valid) can depend on the choice of the marked particles i and j [3, 20, 27]. For the classes 
of models considered in this paper we will be always marking only massless particles, and 
in such a way that the helicities of the marked particles can take the values, 

(h^hj) = (+,-) , (+,+) , (-,-) (2.8) 

but not (hi, hj) = ( — , +). In addition to this, if the two marked particles are a quark and 
an anti-quark of the same flavour, they should not be adjacent. Finally, for adjacent quarks 
and gluons of equal helicity cases we make the choices (jq,ig~) or (iq,jg)- 

3. Massive Vector Bosons 

The principal goal of this paper is to apply on-shell recursion relations to tree-level am- 
plitudes involving massive or off-shell particles with spin coupled to massless particles. To 
begin with we will consider massive vector bosons coupled to massless gauge fields and 
fermions. Essentially we consider a generic theory with a non-Abelian gauge group being a 
product G\ x G2, where G\ is unbroken, and G2 is broken by the Higgs mechanism. Gauge 
fields of the G\ group are massless and the gauge fields of broken group G2 are massive or 
off-shell vector bosons V. Two gauge groups are coupled to each other via fermions which 
are charged under both groups. We will use the 'colour decomposition' representation for 
scattering amplitudes with respect to both groups, and hence the colour-stripped ampli- 
tudes will be purely kinematic quantities which will not depend on the choice of G\ and 
G2 nor on the representations for the matter fermions. 

This set-up is rather general, and in particular it incorporates the elements of the 
Standard Model. In this case G\ is SU (3) and the corresponding gauge fields are gluons 
g; the gauge fields of the (partially) broken group, G2 = SU(2) x 17(1), are massive or 
off-shell vector bosons W ± , Z° and 7*. The fermions can be taken to be (anti)-quarks, 
q, q, transforming in the (anti)-fundamental representations of both groups. Even in the 
general case, we will continue denoting massless gauge fields as gluons, and fermions as 
quarks. Massive vector bosons will be denoted as Vs. 

The quantities we want to consider are the G\- and G2-colour-stripped purely kinematic 
tree-level amplitudes 

S lll ... flm (l q ,2,^,... ,n-l,n q ). (3.1) 

These are the (m + n)-point amplitudes with m massive vector bosons , . . . V^ m coupled 
to n massless partons. More specifically, we consider the case of a single quark-antiquark 
pair 2 denoted in (3.1) as l q , nq and n — 2 gluons labelled 2, 3, . . . , n — 2. 

2 One cannot have less than one qq-p&xr in amplitudes coupling V's to g's at tree level. Amplitudes with 
more than one qq-pair will not be considered in this paper, but they can be calculated in a similar way. 
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The group-theoretical dependence in the amplitudes can be easily restored in the 
usual way. For the case of fundamental fermions the amplitude (3.1) is multiplied by 
(T a2 . ■■T an - 1 ) ilin and by (T bl . . .T bm ) kqk -. Then it is summed over all permutations of 
02, • • • , a n _i and over all permutations of bi, . . . , b m . Here T a and T b are the generators of 
the G\ and the G2 groups respectively. 

The physical states corresponding to all massless particles in amplitudes (3.1) will 
always be represented in the helicity basis, e.g. 1~ , 2 + , 3~ . . . , (n — 1) + , rig- . At the same 
time, for the massive (off-shell) vectors , • • • V^ m we will always choose to not multiply 
them by external wave functions, and instead of helicities or polarisations they will be 
characterised by their Lorentz indices //i,...,// m . Thus, the amplitudes S^...^ are the 
multi- vector boson currents. 

Working with multi-currents (3.1) will first of all facilitate our calculation: single vector 
currents will be used in calculations of double currents and so on, as will be seen in section 5. 
Furthermore, multi-currents can be easily used to calculate full physical amplitudes which 
include the decay of the massive (off-shell) vector bosons into light stable states. This is 
achieved by contracting each Lorentz index \x in (3.1) with the current describing the 
relevant decay mode of each vector boson V^. 

In the Standard Model, for example, one can consider decays of unstable vector bosons 
into a fermion-antifermion (lepton or quark) pair, so that for virtual photon 7* or for 
V = W ± , Z boson decay we have: 

L T - -e Q f Qq -2 

7* 



V V;H V V;H U (Pfh^ U (Pf) 



L V = ~ e ^2 , ^ ,T N • ( 3 - 2 ) 



Here the couplings vy-H for V either W or Z bosons with either left (L) or right (R) handed 
polarisations are given by 



.,// .n^ // _ 0f I l,r^ 2 ^Qf »ih _ 1_ 

J Z;R — ^ifoosOw ' U Z;L— V/ sin^cos^ ' U W;L ~ ^2 sin 



(3.3) 

and all others zero. Uij is the CKM mixing matrix and the rest of notation is standard. 



V W;L ~ ^ sin e w ij ' V W;L ~ ^2 shldw U ji ' V W;L ~ ^2 sin w U n " 



4. Single Vector Boson Currents 

The single vector boson currents were previously calculated by Berends, Giele and Kuijf 
[29] using the recursive technique based on iterations of classical equations of motion [28] . 
More recently these single vector currents were also discussed and derived in [15] using a 
combination of Berends-Giele recursion relations and the MHV rules of [1]. Here we will 
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employ the BCFW on-shell recursion relations of section 2 to derive slightly more compact 
expressions as well as new results for n-parton single currents for some specific helicity 
arrangements of partons. As mentioned earlier, single off-shell currents are not only 
interesting on their own right, more significantly, they play an important part in a recursive 
construction of currents with two and more vector bosons as will be explained in the next 
section. We note that a similar observation has also been made in [15]. 

We now proceed to construct the single vector boson currents, 

S^l x q ,2 h \...,(n- l) h -\n- q x ). (4.1) 

The on-shell recursion relations construct amplitudes from on-shell amplitudes with fewer 
particles. Assuming that one can always avoid contributions at z — ► oo and using the 
recursion relation n — 3 times gives a representation of the n-point amplitude entirely 
in terms of the 3-point vertices. Hence, 3-point vertices are the building blocks of all 
larger amplitudes and hence the on-shell recursion relations reduce the task of computing 
general amplitudes to the computation of all 3-point vertices. This is indeed the case in 
theories with massless vectors coupled to fermions and also to massless and massive scalars 3 
[2, 3, 27]. 

The three-gluon and the quark-gluon-antiquark primitive amplitudes are given by the 
standard MHV and MHV expressions: 

^(l-,2-,3 + ) = , ^ 3 (l + ,2+,3-) = -Jj|EL , (4.2) 

^ 3 (l-,2-,3+) =<gf, ^ s (l7 > 2+,3+) =-^, (4-3) 

.4,(1+ 2" 3^) = , .4,(1+ 2+, 3.) = . (4.4) 

Here as always, all momenta ki are assumed to be complex which ensures that these 3-point 
amplitudes do not vanish on-shell [2, 3]. 

In addition we need to introduce two new primitive vertices for an off-shell vector 
boson coupled to a qq pair. They are derived directly from the Feynman rules: 

^(1-, 2+) = (1|^|2] = [2|<7„|1> , 

^(1+23) = [1|<7„|2> = <2|<7„|1] . (4.5) 

Here, a^aa and a^ a are the standard four Pauli matrices. Our conventions for spinor 
contractions are summarised in the Appendix. 



3 This is however not always the case in theories involving scalar self-interactions. The 4-point vertex 
corresponding to a cj> 4 interaction clearly cannot be reduced to 3-point vertices. Recursion relation cannot 
be applied to this vertex since it is a constant. Hence the corresponding amplitude Ai(z) does not depend 
on z, and this necessarily leads to a non-vanishing contribution at z — * oo. We thank George Georgiou for 
pointing this out to us. 
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2 l q 2 

+ \ 3 

3 » " 4 ? - 

Figure 1: Decomposition of the five-parton amplitude after applying the recursion relation. 

Action of parity symmetry reduces the number of independent currents. Parity trans- 
formation is simply the complex conjugation in terms of S a p := S^a^, 

S ap (l x q , 2 h *,...,(n- l) h ^ , n~ q x ) = (Sp & (l-\ 2~ h * , . . . , (n - l)~ h ^ , n})) * . (4.6) 
This formula is generalised to multi-vector boson currents in an obvious way: 

SaiPi...a m m = (Sp 1 a 1 .../3 m a m )*- ( 4 - 7 ) 

4.1 Single Currents with n = 4, 5,6 partons 

In the four-parton case the recursion relation for reduces to, 

S li (l-,2 ± ,3+) = S li (l-,-^)—A(Pg,2 ± ,3±). (4.8) 

We now have to choose which of the marked particles, 1, 2, is i and which is j. In order to 
avoid boundary contributions according to (2.8) we must choose i = 2 for the 2 + helicity 
and j = 2 for 2~ helicity. This results in, 

W' 2t - S ^-W' (49) 

W , 2 -,3^IM, (4,0, 

where Py is the momentum of the vector boson. Momentum conservation implies Py = 
—Pi — P2 — P3 ■ The two remaining helicity configurations can be obtained by parity trans- 
formations. 

Figure 1 shows the decomposition of the five-parton current. Just as in the four-parton 
case here we mark the quark 1 and adjacent gluon 2, 
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g„(l-2+,3+,4+)= ,\ZoZZ > 



5^(1- 2-, 3+, 4+) 



5 M (l-,2+,3-,4+) 



(12) (23) (34) ' 
[13][3|(l + 2)iV^(l + 2)|3] 
si 23 [12][23](4|2 + 3|l] 
K3 

[24] 3 (l|<r„iV|l) 
s 23 4[23][34](l|2 + 3|4] 



(24)(2|(3 + 4)W(3 + 4)|2) 

5234(23)(34)(4|2 + 3|1] ' 1 ' ) 



(13) 3 [4|JV^|4] 
s 123 (12)(23)(l|2 + 3|4] ' 1 • j 

The S' / i(+; ++; — ) configuration can be obtained from eq. (4.11) by using the "line reversal" 
symmetry, 

2 h \...,n- , n^) = {-l) n +%(n-\ n - l"™" 1 2 h \l x ) . (4.14) 

The other 4 helicity amplitudes can then by obtained via parity transformations. Notice 
that if we use both (4.6) and (4.14) we can also relate (4.10) and (4.9). 

The 6-point amplitudes can be computed in much the same way. We choose to mark 
massless particles in such a way as to generate the most compact analytic expression. The 
following amplitude was computed using i = 2 and j = 1: 

W .2^.3 + .4 + ,5t) = - (i ^W 5> . (4.5) 

The following expressions was derived using i = 3 and j = 4: 

<25><2|(3 + 4 + 5)P v a ll (3 + 4 + 5)|2) 



5 /i (l-,2-,3+,4+,5+) = - 
+ 
+ 



q ' S23 45<23)(34)(45)(5|Py|l] 
[13][3|(l + 2)<yfV(l + 2)|3] 
s 123 (45)[12][23](4|2 + 3|l] 
(2|3 + 4|1](2|(3 + 4)(1 + 2 + 3 + A)a ll P v {l + 2 + 3 + 4)(3 + 4)|2) 
si 234 S234(23)(34)(5|Py|l](4|2 + 3|l] 

(4.16) 



For the choice i = 2 and j = 1 we find: 

(13) 3 (35)(3|(4 + 5)Pya M (4 + 5)|3) 



5 M (l-,2+,3-,4+,5+) = 
+ 



(12)(23)(34)(45)(5|Py (1 + 2)|3)<3|(4 + 5)Py|l) 
(35)(3|4 + 5|2] 3 (l|q M P v |l) 
S2345S34 5 (34) (45) (3 1 (4 + 5)P V |1)(5|3 + 4|2] 
(13) 3 (3|1 + 2|4][4|(1 + 2 + 3)^Py(l + 2 + 3)|4] 

si 23 4Si 23 (12)(23)(5|Py(l + 2)|3)(1|2 + 3|4] 

[24] 4 (l|^Py|l) 



s 23 4(5|3 + 4|2](l|2 + 3|4][23][34] ' 



(4.17) 
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Finally using i = 3 and j = 2 we find, 



S „(l-2+, 3 +,4-,5t>- <">'MV',W 



Sl23 4(12}(23}(34}(l|Py|5] 

[35] 3 (l|a M Py|l) 
S34 5 (12}(2|3 + 4|5][34][45] 



(4|2 + 3|5] 3 (l|^Pyll) 
S234S2345(23)(34)(2|3+4|5](l|Py|5] " l " j 



Using the parity and line reversal symmetries of eqs. (4.6) and (4.14) we can easily obtain 
expressions for the other 12 helicity configurations. 

All the amplitudes presented in this section have been numerically checked against 
Feynman-diagram based calculations. 



4.2 n-point Currents 

It is also possible to construct single vector boson currents with n partons in the helicity 
configurations with maximal helicity violation, next-to-maximal helicity violation and be- 
yond. The current for vector boson decaying to a quark pair and any number of positive 
helicity gluons has been known for some time [29], 

^(1-2+, . . . , (n - 1) + , n+) = (-If niyi^ • ( 4 - 19 ) 

As usual Py is the momentum of the vector boson and Py = —p\ — . . . — p n . This can be 
easily proved by induction using on-shell recursion relations. The fact that any pure QCD 
amplitude with less than two negative helicities is zero guarantees that the only contribution 
to the n-point current involves an (n — l)-point current and an on-shell (complex) 3-gluon 
vertex, as shown in figure 2. This is the first non-vanishing helicity amplitude and hence 
it is the MHV current. For completeness, the other MHV-type currents are given by, 
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\ + - 



1+ 



J=3 



n-l 

+ n * 



n 



1+ 



+ - 




Figure 3: Decomposition of the n-point vector boson current with n — 3 positive helicity gluons 
and two adjacent negative helicities. 



2+ ...,(n-l)+ng 



Sfi(lq ,2 , ... ,(n - 1) ,n q ) = 
^(l-,2-,...,(n-l)-,n±) = 



IE=i<a« + l>' 
[n\P v a^\n] 



(4.20) 



It is interesting to note that eq. (4.19) allows us to immediately write down compact 
expressions for the NMHV currents with both adjacent and non- adjacent minuses. If we 
mark the two negative helicity particles then each sub-amplitude in the recursion relation 
will contain at most 2 negative helicities. Figure 3 shows the decomposition into a sum 
of sub-diagrams. We draw a pure QCD amplitude on the right of each diagram. It must 
contain at least 2 negative helicity particles and this fixes the helicity on the right of the 
propagator to be negative. Helicity conservation then ensures that the vector current on 
the right has only one negative helicity, the marked quark, and so is an MHV current. We 
can therefore use (4.19) to write down the NMHV current, marking i = 1 and j = 2: 



,S M (l-,2-,3 + ,...,(n-l) + ,4) 



-(-1)" / (2|K 3 , n Py<vFC 3 ,n|2}<2n) 



K=2< 



1 [ ■•■ ! '• ' '• 



+ 1)1 s 2 ,n{n\K 2 , n -i\i\ 



n-l 

+ £ 

3=3 



(2\K 3 jK 1 jP v a^K 1 jK 3 j2){2\K 3J \l}{jj + 1) 



S2jSlj01^2J-l|l]0' + l|^2,j|l] y 

This expression is the n-parton generalisation of equations (4.12) and (4.16). 



(4.21) 



We can also consider the case where the helicity along the quark line is flipped. This 
is a special case as we can still eliminate all contributions from NMHV vertices. The result 
is, 



2+,..., (n- 2)+, (n-l) 



■1) 



( (In - l) 3 [n\a^P v \n] 
X[T = \(aa + l)\ si, n -i<l|iVM 



n— 3 



(n - l\K ]+l ^ 2 \nf(l\<j^P v \l){j j + 1) 



J=1 Sj+l,nSj+l,n-l(j + l\K j+ltn -i\n](j\K j+ltn -i\n] J ' 
matching equations (4.13) and (4.18) when n = 4 and 5 respectively. 



(4.22) 
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Figure 4: Decomposition of the n-point vector boson current with n — 3 positive hclicity gluons and 
two negative helicities, separated by a single gluon with positive hclicity. The second contribution 
involves the NMHV current with adjacent negative helicities given in (4.21). 



For NMHV currents with non-adjacent negative helicities we can re-use the above 
result to find the amplitude where the negative helicities are separated by one positive 
helicity. The corresponding diagrams are shown in Fig. 4 where we mark % = 2 and j = 1, 

S „(i-,2+,3-,4+ ...,„-r>±) = _ _ .J^'t^^M 



S2,nS3,n(M K 2,nK4:,n\3)(n\K 3 ^i\2} Y\a=3{ a a + 



n-1 



(l|^P y |l)(3|K4, J |2] 4 (jj + l) 



+ £ 8 2 js 3 j{l\K 2 jK 4J m\K 3 j- 1 \2](j + 1\K 3J \2] m=3(«« + 1) 



Substituting eq. (4.21) and simplifying the shifts results in the following expression: 

1 / (l|Py^|l)(3|K 4 ,n|2] 3 (3n) 



^(l-,2+,3-,4 + ,...,n-l+,n±)- — J 

lla=3 

^4,,|2] 4 (jj 

S2^ 3 {l\K 2tl K^\?,){j\K 3 ^ 1 \2\{j + l\K 3d \2] 



Ul= 3 {aa + 1) V S2,nS3,n(l|^2 ! n^4,n|3}(n| J Fr 3 ,n-l|2] 



(l\<j^P v \l)(3\K^\2f(jj + l) 



n-1 

(l3f{n3)(3\K^ n P v a^ n \3) 



(l\K 2tn K^ n \3)(n\K 1 ^ 1 K 12 \3)(12)(23) 
^ (13)^3\K 4 jK 1 jP v ^K 1J K^\3)(3\K 12 K A j\3)Uj + ^ \ 

s^m^K^m^K^ij + i|^ 1J ^ 12 |3)(i2)(23) y 1 • ) 

Remaining NMHV currents can be constructed in a similar way. We can keep adding an 
extra positive helicity separating the negative helicities. We have checked that the n-parton 
result in eq. (4.24) agrees with eqs. (4.13) and (4.17) for n = 4 and n = 5 respectively. 
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j + l+ 



Figure 5: Decomposition of the single vector boson current with three consecutive negative he- 
licites. The recursion re-uses the current with two consecutive nchativc helicites. 



As a final example we consider the NNMHV current with three adjacent negative 
helicities. 

By marking particles i = 2 and j = 3 we ensure that only the NMHV current (4.21) is 
needed. Explicit calculation yields, 



^(1- 2-,3-,4+,...,n-l+,^ 



(-If 



m=3<««+i> 



(3n)(3jiY 4 , w Py^4,n|3) 

S3 ,n[12](n|K 3 , n -i|2] 



n-1 

+ £ 

j=4 



■i/v;..; ^ •'' K ;.j K l\ (T ii k \.j ^ \.j •"' ././ + 1) 
S3 j SlJ [12]0- + l|^ 2j |l](i|if3j-i|2] 



n— 1 n—l 



si^ 2 ^ 2jS3 j(Z|^-i|l](/ + l|^2,/|l](j + l|K 3 ,i|2](j|K 3 ^i|2] 



EE 

i=4 z=i+i 

" _1 (3\K 4 jK 2j \n)(3\K AiJ K 2 jK 2:n P v a^K 2:n K 2>J K 4:J \3)(jj + 1)' 



E 

J=4 



SZ,jS2,nS2,M K 2,n-l\l](j + l|^ 3> j|2]0'l^3j-l|2] 



(4.25) 



We have explicitly checked eq. (4.25) for up to six partons. 

By repeated use of the recursion formulae, further n-point currents may be obtained. 



5. Double Vector Boson Currents 



We now turn to double vector boson currents S^ v . We start by considering the small- 
est amplitude of this type, the one with only two partons, S^q^q). One might expect 
that on-shell recursion relations can be used to derive S liu (q,q) from two single vector 
boson amplitudes S^q,^). However there is a difficulty in writing down such a recursion 
relation. We cannot mark the two massless particles in S^ u (q,q) since it is known that 
marking adjacent massless quarks results in a non- vanishing boundary contribution [20] to 
the amplitude. Choosing to mark massive particles also leads to (unnecessary) technical 
complications. It is actually much simpler to derive the four-point amplitude S^ u (q, q) from 
Feynman diagrams and use this four-point amplitude as a new primitive vertex in further 
recursive calculations of S^^g, . . . , rig). 
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Figure 6: The two Feynman diagrams contributing to the amplitude with two massive vector 
bosons and a quark pair. 

In fact, we will use a more elegant approach. In general, there are two Feynman 
diagrams contributing to S^ u (q, q), as shown in figure 6. We could evaluate both diagrams 
and use the whole amplitude as a building block for larger amplitudes, however, it is much 
more efficient to split the calculation into two parts in order to re-use the single vector 
boson currents computed in section 4. 

The first part corresponds to the second diagram in figure 6, it contains the non- 
Abelian three-vertex of massive vector bosons. We can compute such contributions to a 
generic S^ v {\ q , . . . ,riq) by contracting two single vector boson currents . . . , n'q) of 

the previous section with the colour ordered Feynman three-point vertex. This approach 
was used to calculate the non-Abelian contribution to S^ u (q,g,q) in reference [15]. Note 
that if one is dealing with uncharged gauge bosons which have no self-coupling, for example 
Z bosons, this term is trivially zero. 

The second part does not contain a non-Abelian three-vertex of vector bosons, it 
corresponds to the first diagram in figure 6. This second Abelian contribution to a generic 
Sfj, v {\q, . . . ,riq) can then be evaluated using on-shell recursion relations. 

Thus, guided by figure 6 we represent the colour ordered double current with n partons 
in the form: 

V(l ? , ...,nq)= T$ p (P Vl , P V2 , -P) (p2 ^ M 2 } S p (l q , • • • , Uq) 

+ S^ dian (l q ,...,nq). (5.1) 

Here 

(5.2) 

is the colour-ordered three-vertex of massive vector bosons, with all momenta defined to 
be in-going. 

The primitive vertex for the Abelian contribution is given by the first Feynman diagram 
in figure 6 which evaluates to: 

S MeUan {l - 2±) = _J_ {1 |^ (1 + p^) | 2] . (5.3) 
SlP V2 
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Figure 7: Recursive decomposition of the five-point Abelian amplitude for two vector bosons, a 
quark pair and a gluon using the recursion relation. 



The remaining non- Abelian part of this four-point amplitude is determined by the first 
line of (5.1) 



^non— Abelian 



(l q ,2 q ) = T$ p {P Vl ,P V2 ,-P) — S<>(l q ,2 q ) (5.4) 



where S p (l q , 2 q ) is given in (4.5). 

In general one needs to determine only the Abelian components of the n-parton double 
currents S^ u (l q , . . . , n q ), the non-Abelian components are fully determined by the first line 
of (5.1) in terms of the known single currents. 

Abelian components are characterised by having massive vector bosons only on external 
lines, and they can always be calculated recursively. In general, in order to calculate the 
Abelian part of any double current, Sffl eban (l q , 2, . . . , nj), one needs to draw all recursive 
decompositions of this current such that the internal line is a quark or a gluon and not a 
massive vector boson. 

First we calculate the two five-point amplitudes by marking the quark and adjacent 
gluon. The recursion relation for S$ elian (l q , 2, 3 5 ) is depicted in figure 7. It yields, 

[13][3|^(l + Py 2 )^(l + 2)|3] 



^Abelian 



(1- 2",3+) = + 



+ 



[12][23][3|P yi Py 2 |l] 
(2|<7 / ,|3][l|g„(l + 2)|3] 

[12][3|P yi Py 2 |l] 
<2|(l + JV 2 )a I/ iV 2 (l + iV 2 )|2)[3|iV 1 g M |3] 

siP V2 S3P Vl [^\Pv 1 Pv 2 m 
(13)<1|(2 + 3)^(1 + Py 2 )a i ,|l) 



(5.5) 



o Abelian t-\ — 9+ o+\ 

""" 1 9 ' ' ~ q) (12)(23)(l|iV 2 iV 1 |3) 



+ 



(1|^|2](1|(2 + 3)^|1) 

(12)(l|Py 2 P yi |3) 

(l|Py 2 ov|l)[2|(l + Pv 2 )P Vl ^(1 + Pv 2 ) 1 2] 



(5.6) 



SlPv 2 S3P Vl (l\Pv 2 Pvi\^) 

The Sfj, u (+; ±; — ) configurations can be obtained from eqs. (4.11) by using either parity 
(4.7) or the "line reversal" symmetry, 



o (-[X nh 2 
> J fJ,U\ J - q ■> ^ > 



,n 



1 



h "-\n7 X ) 



(-lf +1 ^(n- A 



n 



1 



1 nh 2 iA\ 

, . . . , L , lg) 



(5.7) 
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Finally we give results for the six point amplitudes n (l g ,2,3, 4^). Taking the 

generalised parity relation (4.7) and the line-reversal identity (5.7), there are three inde- 
pendent helicity configurations. Again we use on-shell recursion relations and mark the 
quark and adjacent gluon. Choosing i = 2 and j = 1 we find, 



+ 
+ 
+ 



- 2+ 3+ 4+^ = (14)(lk,|2](l|(2 + 3 + 4)q f< |l) 
*' ' ' (12)(13)(34)(l|iV 2 iV 1 |4) 
(14)(l|(2 + 3 + 4)^(l + iV 2 )^|l) 

<12)<23)<34)<l|iV 2 iV 1 |4> 
(l|(2 + 3)^|l)(l|(iV 1 +Py 2 K|l) 

<12)<13)<23)<l|JV 2 iV 1 |4) 
(l|JV 2 g y |l)(l|(2 + 3)(4 + fV 1 )Py ay (4 + P Vl ) (2 + 3) 1 1) 
s 4Pvi (12) (23) (l|iV a P Vl |4> (l|Py 2 (4 + P Vl ) |3) 

(l|Py 2 ^|l)[2|(l+ Py^Py^l + Py,)^] 

^ 2 5l2P v - 2 (34)(l|Py 2 (4 + Py 1 )|3) ' l " j 



and 



■Abelian n - 9+ o- 4+ x (13) 3 (3[1 + 2|4][4|^(1 + Py 2 )^(l + 2 + 3)|4] 
{q, Sl23 (12}(23}(l|2 + 3|4](3|(l + 2)Py 2 P yi |4] 

(13) 3 (3|(4 + P Vl )^P V2 (4 + Py 1 )|3)[4|Py 1 ^|4] 
S4 p Vl (12)(23>(l|Py 2 (4 + Py 1 )|3)(3|(l + 2)Py 2 Pyj4] 

(13) 2 (1|^|4](3|1 + 2|4][2|<r„(l + 2 + 3)|4] 

Sl23 (12)(l|2 + 3|4](3|(l + 2)Py 2 Py 1 |4] 
(13) 2 (3|^|4](3|(1 + 2)^(l + 2 + 3)|4] 

Sl23 (12)(23)(3|(l + 2)Py 2 Py 1 |4] 
[24] 3 (1|3 + 4|2](1|(2 + 3 + 4)<r„(l + Py 2 )a v \l) 
S234(l|iV a iVi(3 + 4)|2](1|2 + 3|4] [23] [34] 
[24] 3 (l|^|2](l|(2 + 3 + 4)^|l) 
(l|Py 2 Py 1 (3 + 4)|2](l|2 + 3|4][23][34] 

[24] 3 (l|Py 2 ^|l)[2|(l + Py 2 )Py 1 q At (l +Py 2 )|2] 

s 1Pv2 (l\P V2 P Vl (3 + 4)|2] [23] [34] [2|(1 + P V2 )P Vl |4] 

(Sll + Py.l^UlPy.^l^KlPy,^!^ 

S4P VlS iP V2 Si2P V2 (l|Py 2 (4 + Py 1 )|3)[2|(l + Py 2 )Py 1 |4]' 1 • j 



- 14 - 



+ 



For the last amplitude we choose i = 1 and j = 2 which yields the following expression, 

q AbeUan n - ^ ^ . + , = (24) [1 1 Q v (P Vl + Py 2 ) (3 + 4) 1 2) (2 1 (3 + 4)<T, |2) 

^ {q ' ' ' ~ q > (23)(34)(2|(3 + 4)P yi P F2 |l](4|2 + 3|l] 

(24)(2|3 + 4|1](2|(3 + 4)^(1 + Pv 2 )a v {P Vl + Py 2 )(3 + 4)|2) 

s 234 (23) (34) (2 1 (3 + 4)P Vl Py 2 | 1] (4 1 2 + 3 1 1] 
(24)(2|(1 + Pv 2 )a v Pv 2 (\ + Py 2 )|2)(2[(3 + 4)P Vl a M (3 + 4)|2) 

s lPv2 (23) (34) (2| (1 + P V2 )P Vl |4) (2| (3 + 4)P Vl Py 2 I 1] 
[13][3|1 + 2|4)[3|(1 + 2)(P yi + Pv 2 )a^l + Py 2 )a v {l + 2)|3] 

si 23 (4|2 + 3|1][12] [23] [3|(1 + 2)Py 2 P yi |4) 
[13] [3 1 ( 1 + 2)P V2 a v ( 1 + 2) 1 3] [3 1 (4 + P Vl )P Vl (4 + P Vl ) 1 3] 
s 4Pvi [12][23][3|(4 + P Vl )Py 2 |l][3|(l + 2)P V2 P Vl |4) 

(2|l + JV 2 |3](2|(l+iVj^JV 2 (l + iVj|2)[3|(4 + iV 1 )iV 1 ^(4 + JV 1 )|3] 
s 4 p Vi s 1Pv2 s 12 p V2 {2\(l + P V2 )P Vl \4)[l\P V2 (4 + P Vl )\3} 

[3|(1 + 2)(P Vl + iV 2 )q M |2)[3|l + 2|4)[1|^(1 + 2)|3] 

s 123 (4|2 + 3|1][12][3|(1 + 2)Py 2 P yi |4) 
[31^(1 + 2)|3][3|(1 + 2){P Vl + Py 2 )a^l + 2)|3] 

Sl23 [12][23][3|(l+2)P y2 P yi |4) • l °- Uj 

The procedure described here can straightforwardly be generalised to processes involv- 
ing three or more vector bosons. In each case, there will be a mixture of terms that either 
involve a triple or quartic gauge boson vertex (non-Abelian) or a new multi-gauge boson 
current (Abelian). The non-abelian contribution is straightfoward and involves currents 
with coupling that couples currents involving fewer gauge bosons. These are in principle 
known. For each additional vector boson, the Abelian contribution must be recomputed. 
There will be a new primitive vertex which can be obtained directly from the single (colour- 
ordered) Feynman diagram. 

This is illustrated in Fig. 8 for three vector bosons. The first diagram yields a new 
primitive vertex S^ lfl2 ^ 3 (l q , 2q) which forms the seed for recursively calculating the Abelian 
contribution to the amplitude. The other three (Non-Abelian) graphs can be straightfor- 
wardly obtained by reusing the single and double vector boson currents. The colour ordered 
triple current with n partons is thus, 

•SVi£t2^3(l<?' • • • ' n <?) = < ~ViM2A t 3 (lg> • • • ' n <?) 

+ T(% 2p (P Vl , Pv 2 , -P12) 2 2 S"« (1„ • • • , ng) 

\ 12 p 12 / 

+ T£IJP V2 , Pv 3 , -P23) (p2 2 ; ■ • • - «ff) 

v 23 P23 / 

+ T^ sp (P Vl ,Pv 2 ,Pv 3 ,-Pi2 3 ) (p2 2 ; S"{l q , ...,n q ). (5.11) 

V 123 P123 / 

The colour ordered quartic gauge boson vertex is given by, 

-^/ii)t2/i3^4(Pi'^'2,P3,P4) — 2g^ il ^ 3 g^ 2 ^ 4 — 1/13 5/i2 M4 — 9hiH29h3H4- 

(5.12) 
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Figure 8: Contributions to the amplitude with three massive vector bosons and a quark anti-quark 
pair. 



6. Recursion Relations for Massive Particles with Spin on Internal Lines 



So far we have been considering application of recursion relations where massive particles 
with spin were absent from the internal lines. In other words, we have been able to set up 
the recursive calculations of double vector boson currents in such a way that the massive 
vector bosons were playing the role of 'external sources' in the left and in the right hand 
vertices, but were not propagating through the recursive diagram. We now would like to 
show how to use recursion relations also for propagating massive particles. In our earlier 
work [27] we have accomplished this for massive scalars, and now we want to generalise 
this approach to massive particles with spin. 

The main difference between internal massive scalars of Ref. [27] and internal massive 
fermions of vector bosons is that the latter have more than one polararisation or spin state. 
In the standard recursion relation (2.1) all particles are assumed to be in a state with fixed 
helicity, and there is a summation over all these states. We want to avoid using helicity 
states for internal massive particles and instead to use a more natural basis of states 

In this section we will describe a new way to implement the recursion relation in 
this case and will illustrate its use by calculating a simple amplitude of two heavy quarks 
scattering into two gluons. 

The main point here is that the sum over helicities h of the internal particle in the 
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standard recursion relation, 

An{Pl, ■ ■ ■ ,Pn) = ^2 ^2A L (Pr, ...,Pi,...,p s , ~P h ) p2 ~ m2 
partitions h P 

X A R {P~ h ,p s +l, ■ ■ ■ ,Pj, ■ ■ ■ ,Pr-l) , (6.1) 

can be replaced by the sum over all of the spin states, rather than helicity quantum numbers 
which are not well suited for massive particles. So, we first replace the sum over helicities by 
the sum over appropriately defined spin states. For massive fermions this is the conventional 
spin sum: 

^2 u s (p)u s (p) =p + m p (6.2) 

s=l,2 

^2 v s (p)v s (p) =f-m v (6.3) 

s=l,2 

The remaining spinors and polarisation vectors of the external massive particles can be left 
unfixed and simplified after squaring the amplitude with the spin sums in the conventional 
way. 

Using this in the recursion relation in which a massive quark propagates between the 
two diagrams we have 

^2-A L (Pr;q, ■ ■ ■ ,Pi, ■ ■ ■ ,Ps, ~Pq)^TZ T^R^q^P^ 1 ' •••>£?■,••■ ,Pr-l;q) 



=-A L (Pr;q, ,Ps, - P *)~B2 2^^ P *^+U ■ ■ ■ ,Pj, ■ ■ ■ ,Pr-l;q) (6.4) 

i Trip 

where P* indicates the external spinor wave-function has been stripped off this amplitude. 
In this way we can use the benefits of using the recursion relations to provide reasonably 
compact formulae for amplitudes with massive particles. 

6.1 Example: Calculation of MO-t, 2, 3, 4 f ) 

We now compute the four point amplitude of a top quark pair scattering to two gluons 
as an example of the method described above. We use on-shell recursion relations and 
mark two massless gluons. This leads to a single recursive diagram with a massive fermion 
propagator. We will show that the contribution of this single diagram precisely matches 
the two Feynman diagrams for this process shown in figure 9. With all particles outgoing 
the recursion relation result is, 

4(14,2,3,%) = {p2 l _ m 2 t) "(Pi)l(P2,&) ^> s (P)u s (P)j #p 3 ,6)v(P4)- (6-5) 

Here P = p\ + p 2 is the momentum on the internal line and £2; £3 are reference spinors 
necessary to specify gluon polarisation vectors e^. We will use the Weyl representation of 
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Figure 9: The two Feynman diagrams contributing to the amplitude ^(lt, 2, 3, 4j). 



the Dirac 7-matrices and polarisation vectors, 



If 



a M 





(p,0 



1 



lOfrl 



r(p£) = — |p>[el 



(6.6) 



^(l t ,2",3 + ,4 f ) 




(P 2 - m 2 )s 2 3 




u(p 4 ). 



First we consider the case where the gluons have opposite helicity, A(lt,2 , 3 + ,4f). It is 
convenient to choose £2 = P3 and £3 = P2 so that, 

( m t P 
P mt 

(6.7) 

We choose the marking prescription % = 3 and j = 2 and this ensures that the shifts on 
the polarisation vectors disappear. It can also be seen that the shift in P is also killed by 
either of the two polarisation vectors and hence we can erase all the hats in equation (6.7). 
This is then exactly equivalent to the first diagram of figure 9. It can be easily shown 
using e~(l, 2) • e + (2, 1) = that, with this particular choice of reference momenta, that the 
remaining second Feynman diagram gives a vanishing contribution and so our recursion 
relation result is in agreement with the Feynman diagrams answer. 

The amplitude with both gluons of negative helicity A(lt, 2 + , 3 + , 4j) is of a non-MHV 
type and it provides another interesting test of the recursion relation, which this time 
requires a little algebra. The recursion relation reads: 

' |3)[2|\ (m t P\ ( |2>[3|\ 

i2] ( 3i ){p mt ) [m ) v{p ^ 

(6.8) 

Again, choosing i = 3 and j = 3 removes the shifts on the propagator and the polarisation 
vector of gluon p%. However in this case all the shifts do not vanish as \2] = |2] — z\l] hence 
we are left with the exact expression for the first Feynman diagram plus an extra term 
coming from the surviving shifts: 

' |3)[3|\ (m t f \ ( |2)[3| N 
|3](3| U mt l|3](2| 



.4(l t ,2+ 3+ 4 f ) = 



(P 2 -m 2 )(23} 2 



u(j>i 



(P 2 -m t 2 )(23) 2 
which simplifies to, 

1 



u(pi) 



v{p A ), (6.9) 



(2|p 4 |3](23> 



|3)(2|p 4 |3][3| 
|3](3|p 4 |3]<2| m t |3](32)[3| 



V(p4,). 



(6.10) 
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If the result from the recursion relation is to match the result of the Feynman calculation 
this expression should be equivalent to the second diagram in figure 9. Making use of the 
Dirac equation one can simplify the Feynman calculation to 

~^2«(Pl) ( p 3 q) «(P4)- (6-H) 

It may not be immediately obvious that the expressions (6.10) and (6.11) are equivalent, 
but they are. Firstly we note that the four component spinor can be written in terms of 
two component spinor s, 




(6.12) 



This allows us to expand out (6.10) and immediately identify that the top row is in the cor- 
rect form. The bottom row can be simplified by using mt |i>4] = f 4 \v^) and by decomposing 
\vi) = a|3) + j3\2) we can re-form the bottom row into the correct form and reconstruct 
(6.11). 

This shows that recursion relations can be used to successfully calculate amplitudes 
with massive particles with spin on internal lines, as expected. 



7. Conclusions 

In this paper we have applied the recursion relations derived by Britto, Cachazo, Feng and 
Witten to a range of processes involving both internal and external, off-shell or massive 
particles. Our analysis includes massive particles with spin and this complements our 
earlier work [27] where we have considered massive scalars coupled to massless partons. As 
a first step, in Sec. 4 we derived compact expressions for single vector boson currents where 
an off-shell gauge boson couples to a quark-antiquark pair and up to four gluons. These 
results both agree with and are more compact than the expressions previously available in 
the literature. New results for particular helicity configurations involving any number of 
partons were also derived. 

Adding more off-shell vector bosons is also straightforward. The most direct approach 
is to divide the amplitude into its Abelian and Non-Abelian components. The Non-Abelian 
contribution is straightforwardly obtained by linking one or more currents with fewer vector 
bosons via triple or quartic vertices. On the other hand, the Abelian contribution is 
amenable to a recursive approach using the single Feynman diagram where all of the off- 
shell vector bosons couple to a fermion line in an ordered way. As an example, Sec. 5 
contains compact expressions for the double vector boson current coupling to a quark- 
antiquark pair and up to two gluons in all helicity configurations. 

Note that propagating particles with mass only show up in the Non-Abelian con- 
tribution and we were able to organise the recursive calculation of the Abelian part of 
multi-vector boson currents without having massive intermediate states. However, often 
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it is necessary to deal with propagating massive particles and the recursion relations can 
handle this situation as well. The key point is to generalise the recursion relation to avoid 
using helicities for internal massive particles. The sum over helicities is then replaced by 
the conventional spin sum. 

Taken together, the generalisations of the recursion relations presented here are capable 
of providing an efficient way of calculating tree-level amplitudes within the Standard Model 
and beyond. 
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A. Spinor Conventions and Useful Identities 

We use the standard a^aa and a^ a matrices to represent a vector p^ in spinor notation 

Spinor indices are raised and lowered with e-symbols, such that 

~< = ( A -2) 
We also have an identity a^ a a^ = 26^6^. 

Null vector p^ is represented in terms of dotted and undotted 2-spinors as 

P 2 a = a"a a . (A.3) 

The spinor products are defined by 

(ab) = a a b a = (a\ a \b) a , (A.4) 
[ab] =aab d = [a\ & \bf. (A.5) 

We note that the spinor summation conventions for dotted spinors in (A.5) follow what is 
usually used in QCD literature and differ by a minus sign from [2, 3, 27, 48]. 

Parity transformation requires complex conjugation which for spinors (of a real vector) 
is defined as 

= \a]a => (ab)* = -[ab]. (A.6) 

The spinor sandwiches are defined in such a way that the summation over adjacent 
indices goes from up to down a a for undotted indices and from down to up „" for dotted 
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indices, precisely as in (A. 4) and (A. 5) 



(i\p\j] = X?Paa~Xf, (A.7) 

(i\PrPs\j] = KPraapf X jf3 , (A.8) 

[i\PrP S \j] = wrp SQ/3 A/- ( A - 9 ) 

For massless momenta we have 

Sij = (ij) [j i] , (i \Pr \j] = (i r) [r j] , (A. 10) 

(i\PrPs\j) = (ir)[rs](sj) , [i\p r p s \j] = [ir](rs)[sj] , (A. 11) 

and so on. 

Equation (A.7) also implies that 

(a\a^\b] = (a\ a a^\bf , [b\a»\a) = [b\^ a \a) a (A.12) 
and the identity (A. 2) gives 

(a\a^\b] = [6|<7», (A.13) 



or more generally 



{a\...a^...\b] = [6|...^...|a>, (A.14) 
(a\...a ll ...\b) = -(b\...a^..\a), (A.15) 
[a|...^...|6] = -[6|...^...|a]. (A.16) 



References 



[1] F. Cachazo, P. Svrcek and E. Witten, MHV vertices and tree amplitudes in gauge 
theory, JEEP 09 (2004) 006 [hep-th/0403047]. 

[2] R. Britto, F. Cachazo and B. Feng, New recursion relations for tree amplitudes of 
gluons, Nucl. Phys. B715 (2005) 499-522 [hep-th/0412308]. 

[3] R. Britto, F. Cachazo, B. Feng and E. Witten, Direct proof of tree-level recursion 
relation in Yang-Mills theory, Phys. Rev. Lett. 94 (2005) 181602 [hep-th/0501052]. 

[4] Z. Bern, L. J. Dixon, D. C. Dunbar and D. A. Kosower, One loop n point gauge 
theory amplitudes, unitarity and collinear limits, Nucl. Phys. B425 (1994) 217-260 
[hep-ph/9403226]. 

[5] Z. Bern, L. J. Dixon, D. C. Dunbar and D. A. Kosower, Fusing gauge theory tree 
amplitudes into loop amplitudes, Nucl. Phys. B435 (1995) 59-101 [hep-ph/9409265]. 

[6] C.-J. Zhu, The googly amplitudes in gauge theory, JEEP 04 (2004) 032 
[hep-th/0403115]. 



- 21 - 



[7] G. Georgiou and V. V. Khoze, Tree amplitudes in gauge theory as scalar MEV 
diagrams, JEEP 05 (2004) 070 [hep-th/0404072]. 

[8] D. A. Kosower, Next-to-maximal helicity violating amplitudes in gauge theory, Phys. 
Rev. D71 (2005) 045007 [hep-th/0406175]. 

[9] J.-B. Wu and C.-J. Zhu, MEV vertices and scattering amplitudes in gauge theory, 
JEEP 07 (2004) 032 [hep-th/0406085]. 

[10] J.-B. Wu and C.-J. Zhu, MEV vertices and fermionic scattering amplitudes in gauge 
theory with quarks and gluinos, JEEP 09 (2004) 063 [hep-th/0406146]. 

[11] G. Georgiou, E. W. N. Glover and V. V. Khoze, Non-MEV tree amplitudes in gauge 
theory, JEEP 07 (2004) 048 [hep-th/0407027]. 

[12] V. V. Khoze, Gauge theory amplitudes, scalar graphs and twistor space, 
hep-th/0408233. 

[13] L. J. Dixon, E. W. N. Glover and V. V. Khoze, MEV rules for Eiggs plus 
multi-gluon amplitudes, JEEP 12 (2004) 015 [hep-th/0411092]. 

[14] S. D. Badger, E. W. N. Glover and V. V. Khoze, MEV rules for Eiggs plus 
multi-parton amplitudes, JEEP 03 (2005) 023 [hep-th/0412275]. 

[15] Z. Bern, D. Forde, D. A. Kosower and P. Mastrolia, Twistor-inspired construction of 
electroweak vector boson currents, hep-ph/0412167. 

[16] T. G. Birthwright, E. W. N. Glover, V. V. Khoze and P. Marquard, Multi-gluon 
collinear limits from MEV diagrams, hep-ph/0503063. 

[17] T. G. Birthwright, E. W. N. Glover, V. V. Khoze and P. Marquard, Collinear limits 
in QCD from MEV rules, hep-ph/0505219. 

[18] R. Roiban, M. Spradlin and A. Volovich, Dissolving N = 4 loop amplitudes into 
QCD tree amplitudes, Phys. Rev. Lett. 94 (2005) 102002 [ hep-th/0412265]. 

[19] M.-x. Luo and C.-k. Wen, Recursion relations for tree amplitudes in super gauge 
theories, JEEP 03 (2005) 004 [hep-th/0501121]. 

[20] M.-x. Luo and C.-k. Wen, Compact formulas for all tree amplitudes of six partons, 
Phys. Rev. D71 (2005) 091501 [hep-th/0502009]. 

[21] R. Britto, B. Feng, R. Roiban, M. Spradlin and A. Volovich, All split helicity 
tree-level gluon amplitudes, Phys. Rev. D 71, 105017 (2005) [hep-th/0503198]. 

[22] J. Bedford, A. Brandhuber, B. Spence and G. Travaglini, A recursion relation for 
gravity amplitudes, hep-th/0502146. 

[23] F. Cachazo and P. Svrcek, Tree level recursion relations in general relativity, 
hep-th/0502160. 



- 22 - 



[24] Z. Bern, L. J. Dixon and D. A. Kosower, On-shell recurrence relations for one-loop 
QCD amplitudes, hep-th/0501240. 

[25] Z. Bern, L. J. Dixon and D. A. Kosower, The last of the finite loop amplitudes in 
QCD, hep-ph/0505055. 

[26] Z. Bern, L. J. Dixon and D. A. Kosower, Bootstrapping multi-parton loop amplitudes 
in QCD, hep-ph/0507005. 

[27] S. D. Badger, E. W. N. Glover, V. V. Khoze and P. Svrcek, Recursion relations for 
gauge theory amplitudes with massive particles, hep-th/0504159. 

[28] F. A. Berends and W. T. Giele, Recursive calculations for processes with n gluons, 
Nucl. Phys. B306 (1988) 759. 

[29] F. A. Berends, W. T. Giele and H. Kuijf, Exact expressions for processes involving a 
vector boson and up to five partons, Nucl. Phys. B321 (1989) 39. 

[30] F. Cachazo, P. Svrcek and E. Witten, Twistor space structure of one-loop amplitudes 
in gauge theory, J 'HEP 10 (2004) 074 [hep-th/0406177]. 

[31] A. Brandhuber, B. Spence and G. Travaglini, One-loop gauge theory amplitudes in N 
= 4 super yang-mills from MHV vertices, Nucl. Phys. B706 (2005) 150-180 
[hep-th/0407214]. 

[32] F. Cachazo, P. Svrcek and E. Witten, Gauge theory amplitudes in twistor space and 
holomorphic anomaly, JEEP 10 (2004) 077 [hep-th/0409245]. 

[33] I. Bena, Z. Bern, D. A. Kosower and R. Roiban, Loops in twistor space, Phys. Rev. 
D71 (2005) 106010 [hep-th/0410054]. 

[34] C. Quigley and M. Rozali, One-loop MHV amplitudes in supersymmetric gauge 
theories, J 'HEP 01 (2005) 053 [hep-th/0410278]. 

[35] J. Bedford, A. Brandhuber, B. Spence and G. Travaglini, A twistor approach to 
one-loop amplitudes in N = 1 supersymmetric yang-mills theory, Nucl. Phys. B706 
(2005) 100-126 [hep-th/0410280]. 

[36] J. Bedford, A. Brandhuber, B. Spence and G. Travaglini, N on- supersymmetric loop 
amplitudes and MHV vertices, Nucl. Phys. B712 (2005) 59-85 [hep-th/0412108]. 

[37] F. Cachazo, Holomorphic anomaly of unitarity cuts and one-loop gauge theory 
amplitudes, hep-th/0410077. 

[38] R. Britto, F. Cachazo and B. Feng, Computing one-loop amplitudes from the 
holomorphic anomaly of unitarity cuts, Phys. Rev. D71 (2005) 025012 
[hep-th/0410179]. 



- 23 - 



[39] Z. Bern, V. Del Duca, L. J. Dixon and D. A. Kosower, All 

non-maximally-helicity-violating one-loop seven-gluon amplitudes in N = 4 
super-yang-mills theory, Phys. Rev. D71 (2005) 045006 [hep-th/0410224]. 

[40] S. J. Bidder, N. E. J. Bjerrum-Bohr, L. J. Dixon and D. C. Dunbar, N = 1 

supersymmetric one-loop amplitudes and the holomorphic anomaly of unitarity cuts, 
Phys. Lett. B606 (2005) 189-201 [hep-th/0410296]. 

[41] R. Britto, F. Cachazo and B. Feng, Generalized unitarity and one-loop amplitudes in 
N = 4 super-yang-mills, hep-th/0412103. 

[42] Z. Bern, L. J. Dixon and D. A. Kosower, All next-to-maximally helicity-violating 
one-loop gluon amplitudes in N = 4 super-yang-mills theory, hep-th/0412210. 

[43] S. J. Bidder, N. E. J. Bjerrum-Bohr, D. C. Dunbar and W. B. Perkins, One-loop 
gluon scattering amplitudes in theories with N < 4 super symmetries, Phys. Lett. 
B612 (2005) 75-88 [hep-th/0502028]. 

[44] R. Britto, E. Buchbinder, F. Cachazo and B. Feng, One-loop amplitudes of gluons in 
SQCD, hep-ph/0503132. 

[45] A. Brandhuber, S. McNamara, B. Spence and G. Travaglini, Loop amplitudes in pure 
Yang-Mills from generalised unitarity, hep-th/0506068. 

[46] E. I. Buchbinder and F. Cachazo, Two-loop amplitudes of gluons and octa-cuts in N 
= 4 super Yang-Mills, hep-th/0506126. 

[47] Z. Bern, N. E. J. Bjerrum-Bohr, D. C. Dunbar and H. Ita, Recursive calculation of 
one-loop QCD integral coefficients, hep-ph/0507019. 

[48] E. Witten, Perturbative gauge theory as a string theory in twistor space, Commun. 
Math. Phys. 252 (2004) 189-258 [hep-th/0312171]. 

[49] F. Cachazo and P. Svrcek, Lectures on twistor strings and perturbative Yang-Mills 
theory, hep-th/0504194. 

[50] V. P. Nair, Noncommutative mechanics, Landau levels, twistors and Yang-Mills 
amplitudes, hep-th/0506120. 



- 24 - 



